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Abstract 

Time series of string tension of a simple pendulum has not yet been 
a interesting motion information, even nowadays using a force tension 
sensor can be measured easily. A numerical procedure is presented 
how to obtain motion of pendulum bob from the string tension. Un- 
fortunately, proposed procedure does not work. A quadratic term of to 
loses the sign of u that prevents the procedure to work. Applying sign 
of string tension into the procedure can produce the time series but 
only in positive value. The error for reproduced angular displacement 
is about 7.6 %. 
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1 Introduction 



Pendulum is an interesting system. It seems simple but it does not so. Even 
for a Galileo, it was crucial throughout his career [1]. The simple kind of 
pendulum is known as simple pendulum, which swings in small angle regime 
and has well known periode Tq = 211^1 j g |2l [U H]. Complexity can be 
added to the system of a simple pendulum such as wide angle regime [S], 
elastic string [0], or shaken pivot at where the string attached [7]. Especially 
for wide angle regime, Cristiaan Huygens in 1658 stated that a pendulum, 
which swings in very wide arcs of about 100 °, will introduce inaccuracy 
that is causing the swing period to vary with amplitude changes caused 
by small unavoidable variations in the driving or restoring force [S]. There 
are many attempts to formulate the periode of a simple pendulum for wide 
angle regime, for example, derived from energi consideration [9], based on 
the arithmetic-geometric mean [10], and by introducing a mutliplier [TT] . 
The simple experiment to measure the period automatically is by using a 
photogate sensor which registers only the time as the string is passing the 
sensor [12] . Information of string tension is not used to determine the motion 
but only to confirm the tension calculated from other motion parameters [T3] . 
Based on experiment setup in [H] a calculation procedure to obtained motion 
of pendulum bob is reported in this work. 



2 Theory 



A simple pendulum system consists of a massless string with length / and a 
pendulum bob with mass m as illustrated in Figure [U It is considered in the 
system that air friction has no influence and the string can not be streched. 

Instead of writing the equation of angular position of pendulum bob Q 



^ + fsmO = 0. (1) 

which is derived from Newton's second law of motion in tangential direction, 
the law in radial direction is presented 
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Figure 1: System of simple pendulum with all considered forces: tension 
force F and earth gravitational force mg. 



- +^cos^^- — =0, 2 
at J I ml 

where string tension F and angular position 6 is function of time t. The 
simplest way to solved Equation ([2]) to obtain solution of 6 is by solving 



— = ±\ — ---cosO (3) 
dt \ ml I ^ ^ 

numerically with some initial conditions. Equation ([3]) is a little bit different 
than that is proposed in [9] which is 



f = ±yy(cos^^-coseo), (4) 

with initial conditions ^(0) = +9o and d9/dt{0) = 0. Equation ([3]) can be 
solved numerically. Using Euler method it can be written as 



e±(t + At) = e±(t) ± J ^ - ? cose±(t)At. (5) 

V ml I 

Equation ([5]), for both 9±, will be tested using artificial data produced by 
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Equation ([T]) dan Equation ([T]) will give the following numerical equa- 
tions to produce 9(t) 



a{t) = —— sin6'(t), 

a;(t + At) =uj{t) + a{t), 
e{t + At) = e{t)+io{t) 



(6) 

(7) 
(8) 



with initial conditions 6{t) = 9q and cu = at t = 0. The produced time 
series of 6{t) will be used to produce F{t) using 



3 Results and discussion 

The 9{t) has been simulated with / = 1 m, (yf = 10 m/s^, = 0.5236 rad, 
1^0 = 0, tmax = 10 s, At = 10~^ s, and m = 0.2 kg. Using this time series the 
value of u{t) is calculated using 



and then the string tension F(t) is calculated using Equation (|9]), ( ITOj) . and 
the time series 0{t). 

Figure [2] shows the result of simulated 9(t) and the string force F(t) produced 
from the 6{t) is given in Figure 121 Time series of F{t) has periode double 
than the time series of 9(t). At time t about 0.5 s the pendulum bob passes 
a point where it has minimum height. This point, also at t about 1.5, 2.5, 
and other similar points, divide the range of one oscillation periode into two 
identical region where value of the string tension in these two regions with 
the same \d\ has the same value. That is way the periode of string tension 
F{t) has periode double than the 0{t). In this model only magnitude of F{t) 



F{t) = mlLo'^{t) +mg cos 9{t). 



(9) 



e{t) - e{t - At 

At 



(10) 
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Figure 2: Produced times series 6{t) with I - 
0.5236 rad, Uq = 0, t^ax = 10 s, and At = lO'^ 



1 m, 5( = 10 m/s^, Oq 



s. 



is calculated that the direction always point to a point where the system is 
hung. The mathematical explanation is simpler, Equation (|9]) gives always 
positive value from 9(t): the and cos[6'(t)]. 

Unfortunately, direct application of Equation (|5]) to calculate 6{t) from F{t) 
give unintepreted results which has no physical meaning since value of 6{t) 
can be more than ±2-7? rad. This result is shown in Figure HI This could be 
lied on the ± sign in right side of Equation (|5]). Further investigation shows 
that the calculated 9(t) does not produce the same F(t) when it is used to 
calculate F(t) using Equation ([9]). It means that the integration used in this 
report does not guarantee that the results are self consistent. The different 
between F{t) and F[6{t)] which is represented as AF is given in Figure El A 
correction must be applied to overcome this problem. 

An algorithm that could reduce the value of AF could probably give the 
right 6{t). This is the explanation why the previous result of 6{t) could have 
such a large value 6 G [—6 : 6] rad even the initial angle only 0.5236 rad. 
An iteratif self consisten field scheme is used to overcome this obstacle, but 
unfortunately can not handle this problem. 

Other way is to modify Equation ([5]) by notifying the behaviours of F{t) and 
6{t) from Figure [2] and El Then Equation (E]) is modified into 
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Figure 3: Produced times series F{t) with I = 1 m, g = 10 m/s^, m = 0.2 kg 
and At = 10^^ s from previous 9{t). 



This modification can be implemented since the time series of F{t) is aheady 
available. The result is given in Figure El 

The difference between 6{t) and calculated 6'caic(^)5 which is produced from 
F{t), is defined as 



Ae = \eit)\~9,^Ut). 



(12) 



Plot of A6 against t is given in Figure [71 



4 Conclusion 



A numerical procedure to obtain time series of angular displacement 6{t) 
from string tension time series F{t) of a simple pendulum has been reported. 
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Figure 4: Produced times series 6{t) from F{t) in Figure |3] for: 6^ (o) and 

0- (□). 



The similar Form of the equation that holds for tangential direction must be 
treated differently for radial direction since there is a quadratic term that 
hides the information of direction of angular velocity. A modification to the 
reqular form is shown and proved that it works. Only one minor problem 
remains, that the produced angular displacement has only positive values. 
The difference between orignal and reproduced angular displacement is below 
0.04 rad for initial angular displacement 0.5236, or about 7.6 % error. 
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Figure 5: The difference between F{t) and F[9{t)] for: 6'+ (g) and 9_ (□). 
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Figure 6: Produced 6'caic(^) from F(t) in Figure [31 
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Figure 7: The difference between 9{t) and calculated ^caic(^)- 
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